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Abstract. Hyperfine structure A and B factors of the atomic 5s 5p ®P2 —^ 5s 6s ®Si transition are deter¬ 
mined from collinear laser spectroscopy data of and Nuclear magnetic moments 

and electric quadrupole moments are extracted using reference dipole moments and calculated electric field 
gradients, respectively. The hyperfine structure anomaly for isotopes with Si/2 and ^5/2 nuclear ground 
states and isomeric hii /2 states is evaluated and a linear relationship is observed for all nuclear states 
except S1/2. This corresponds to the Moskowitz-Lombardi rule that was established in the mercury region 
of the nuclear chart but in the case of cadmium the slope is distinctively smaller than for mercury. In total 
four atomic and ionic levels were analyzed and all of them exhibit a similar behaviour. The electric field 
gradient for the atomic 5s 5p ®P2 level is derived from multi-configuration Dirac-Hartree-Fock calculations 
in order to evaluate the spectroscopic nuclear quadrupole moments. The results are consistent with those 
obtained in an ionic transition and based on a similar calculation. 

PACS. 21.10.Ky Electromagnetic moments - 31.30.Gs Hyperfine interactions and isotope effects 


1 Introduction 

Atomic spectroscopy of short-lived isotopes is a unique 
tool to investigate nuclear moments, spins and charge radii 
in a nuclear-model independent way. The magnetic mo¬ 
ments provide access to the single-particle structure of 
nuclei while information on the nuclear shape and col¬ 
lectivity is imprinted in the quadrupole moments. Com¬ 
bined with the spin, extracted from the hyperfine struc¬ 
ture, these nuclear properties and their evolution along 
an isotopic chain provide invaluable insight into the forces 
governing the nucleus. Experimental accuracies of the hy¬ 
perfine structure (hfs) splittings are often at the 10“^ level 
if high-resolution laser spectroscopy is applied. Extraction 
of the nuclear magnetic dipole moment fj, and the spectro¬ 


scopic electric quadrupole moment Qs directly from the 
hyperfine A and B parameters of a particular isotope re¬ 
quires a very good knowledge of the hyperfine fields, i.e. 
the magnetic field Be{0) and the electric field gradient 
Vzz{0), induced by the electrons at the site of the nucleus. 
Alternatively, the moment of the isotope X can be cal¬ 
culated if the respective nuclear moment is known for at 
least one reference isotope of the same isotopic chain. In 
the case of the nuclear magnetic moment, the correspond¬ 
ing relation 

fJ-X = -j-7-MRef, (1) 

TlRef 7 Ref 

can be used to determine the magnetic moment from the 
ratio of the magnetic coupling constants Ax and Afief, 
measured by laser spectroscopy, and the corresponding nu- 




2 


Nadja Frommgen et al.: Collinear laser spectroscopy of atomic cadmium 


clear spins Ix and /Ref - Magnetic dipole moments of stable 
and reasonably long-lived isotopes are usually known with 
high accuracy from nuclear magnetic resonance (NMR) or 
atomic beam magnetic resonance experiments. However, 
Eq. (1) is based on a point-dipole approximation for the 
nuclear magnetic moment. In reality the magnetic mo¬ 
ment is distributed over the finite size of the nucleus. 
The influence of this spatial distribution is a deviation 
of the measured dipole coupling constant A from that of a 
nucleus with identical charge distribution but point-like 
dipole moment, Ap. This is taken into account by the 
Bohr-Weisskopf correction eew according to [1] 


^ • (1 + ebw) • (2 ) 

The correction coefficient eew is isotope dependent and 
reflects the relative contributions of spin and orbital angu¬ 
lar momentum of the active valence nucleons to the mag¬ 
netic moment. As point nuclei do not exist and theoreti¬ 
cal calculations of A factors are not accurate enough, the 
“Bohr-Weisskopf effect” is only seen in small deviations 
of the A-factor ratios between different isotopes from the 
ratios given by the magnetic moments. These deviations 

w Ebw “ ^Bw ( 3 ) 


are known as hyperfine structure anomalies and Eq. (1) 
has to be modified to 


MX 


Ax Ix 

ARef /Ref 


MRef (1 + . 


( 4 ) 


The hfs anomaly is typically of the order of lO”"* for light 
elements to 10“^ for heavy elements and is usually tabu¬ 
lated in percent. 

The study of hfs anomalies requires very accurate and 
independent measurements of magnetic moments and hfs 
parameters. This is the reason why systematic data are 
available only in a few cases. One of the prominent exam¬ 
ples is the series of mercury isotopes between ^^^Hg and 
203 Hg for which Moskowitz and Lombardi [2] realized that 
the hfs anomaly follows a surprisingly simple rule. Plotted 
against the difference of the inverse magnetic moments, 
the anomalies 


1^2 _ ^1 4i M2 
A2 h Ml 


( 5 ) 


reveal a linear dependence corresponding to the relation¬ 
ship 



( 6 ) 


which suggests that for the individual isotopes 


a 

ebw = — ■ 
M/ 


( 7 ) 


Eor the mercury 6s 6p ^Pi level the proportionality con¬ 
stant was found to be a = ±1.0 • 10~^/iAr, depending on 
the coupling I = / ± 1 / 2 , respectively, of orbital and spin 
angular momentum in the shell-model description of the 
odd-neutron state. According to the rule (7) nuclear states 


in the tin region like 1 / 5 / 2 , hii/ 2 , resulting from the coupling 
/ = / ± 1 / 2 , and (i 3/2 from / = / — 1/2 should also exhibit 
equal absolute values of a with opposite sign. 

While a similar behavior was observed for a few iso¬ 
topes of odd-Z elements in the mercury region [3], namely 
I9i,i93jj.^ 196-199^^^ 203,205/p]^^ Conclusive picture about 

the validity of the rule has emerged from scattered data 
on elements in other regions of the nuclear chart [4]. It 
seems desirable to obtain more experimental information 
and to better understand the nuclear model assumptions 
to be made for a theoretical validation of relations (6) and 
(7). Nevertheless, the Moskowitz-Lombardi rule has often 
been used to estimate uncertainties arising from hyperfine 
anomalies for nuclear moments extracted from hyperfine 
splittings (see, e.g., [5,6]). 

The cadmium chain is in many respects similar to the 
mercury case. Cd {Z = 48) is one proton pair below the 
Z = 50 shell closure, as is Hg {Z = 80) below the next shell 
closure at Z = 82. In both regions isomeric states exist 
due to the opposite-parity high-spin orbitals, hn/^ below 
N = 82 and A 3/2 below N = 126. Moreover, accurate mea¬ 
surements of nuclear moments exist for i07,i09,iii,ii3(/<(^ 

[7] , and for ^^®Cd as well as for the isomers 

[8] , covering nuclear spins of J = 5 / 2 , 1/2 and 11 / 2 . This 
facilitates an analysis of hfs anomalies. 

While magnetic moments are usually extracted from 
measured hfs A factors according to Eq. (1) by using ex¬ 
perimental reference values of both quantities for one iso¬ 
tope, a similar approach is not possible for the nuclear 
quadrupole moments since directly measured reference val¬ 
ues of Qs are not available. Instead, the extraction of Qs 
from atomic B factors has to rely on calculations of the 
electric held gradient V^z at the site of the nucleus. Tra¬ 
ditionally such calculations have been based on classical 
hfs theory, using semi-empirical expectation values {r~^) 
for a particular valence electron conhguration. To date, 
all tabulated quadrupole moments of cadmium isotopes 
are taken from the early work of McDermott et ah, in 
particular [9]. The held gradient was calculated from the 
dipole and quadrupole interaction constants of the P 3/2 
electron in the hfs of the ^P 2 and ^Pi states of ^°®Cd, 
and a typical 10% uncertainty was estimated, covering 
the range of values resulting from alternative evaluations 
of Sternheimer-type shielding corrections were ne¬ 

glected. It was a goal of the present work to provide more 
reliable electric held gradients of different atomic states 
from modern hfs theory and thus obtain a reference for 
the evaluation of a new set of Cd quadrupole moments. 

Already in Ref. [10] the quadrupole moments were 
based on a Vzz value taken from Dirac-Hartree-Eock cal¬ 
culations in the Cd+ (Cd H) 5p level. In Sect. 4.5 
we give details of these calculations and of similar ones 
for the atomic (Cd I) 5s 5p ^P 2 level with a comparison 
of the accuracy and consistency of the results. For reasons 
discussed in Sect. 4.6, we rely on the values calculated for 
Cd H and the results given here are consistent with those 
already published in [10]. 

Here we report on laser spectroscopy measurements in 
the atomic 5s 5p ^P 2 —>■ 5s 6s ^Si transition for i06-i24(-;^^ 
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and ^^®Cd. Nuclear ground states as well as long-lived iso¬ 
meric states were investigated. The focus of this paper 
is on the extraction of the nuclear parameters from the 
atomic spectra. Magnetic moments are extracted from hfs 
A factors using Eq. (1) and this is combined with an anal¬ 
ysis of the hfs anomaly and an examination of the appli¬ 
cability of the Moskowitz-Lombardi rule. Quadrupole mo¬ 
ments are extracted from B factors using newly calculated 
electric field gradients obtained from multi-configuration 
Dirac-Hartree-Fock theory. The reliability of these calcu¬ 
lations is verified by comparing the results for the atomic 
level investigated here with those for the level in the Cd’*' 
ion reported in [ 10 ]. 


2 Experimental setup 

The measurements were performed at the on-line isotope 
separator ISOLDE at CERN. The cadmium isotopes were 
produced by 1.4-GeV protons impinging on a uranium 
carbide target and ionized by resonance laser ionization 
in the laser ion source of ISOLDE. After acceleration to 
50 keV the ions were mass separated in the General Pur¬ 
pose Separator (GPS) and guided to the collinear laser 
spectroscopy setup (COLLAPS). At the COLLAPS beam 
line, the fast ion beam is collinearly superimposed with the 
laser beam and subsequently the ions are neutralized by 
charge-exchange reactions with sodium atoms in a vapour 
cell to which the post-acceleration voltage was applied. 
Resonance fluorescence light was observed in the fluores¬ 
cence detection region (FDR), mounted adjacently to the 
exchange cell, by four photomultiplier tubes. The layout 
of the FDR system is described in [11]. 

Spectroscopy was performed in the atomic (Cd I) tran¬ 
sition from the metastable 5s 5p ^P 2 level, partially popu¬ 
lated in the charge-exchange process, to the excited 5s 6 s 
^Si level. This corresponds to a transition wave number of 
19657.031 cm“^ [12]. The required laser light at a wave¬ 
length of 508.7 nm was produced by a dye laser, operated 
with Goumarin 521, which was pumped by an argon-ion 
laser operated at 488 nm. The laser frequency was con¬ 
stantly measured with a wavemeter and locked to an ex¬ 
ternal cavity which was again stabilised in length to a 
frequency stabilised heliummeon laser. Thus a stability of 
the dye laser of the order of a few MHz over several hours 
was obtained. A laser power of about 1 mW was used for 
spectroscopy. 

While the laser was fixed at a frequency in the labo¬ 
ratory system, the velocity of the ions was manipulated by 
applying an additional tunable post-acceleration voltage 
in the range of ilOkeV to the charge-exchange cell. Thus 
the laser frequency is shifted in the rest frame of the ion 
according to the Doppler formula to i^c = i^L 7 (l ~/3). The 
velocity (i = vicin terms of the speed of light is obtained 
from the difference of the (positive) acceleration poten¬ 
tial at the ISOLDE ion source and the post-acceleration 
voltage at the charge-exchange cell. The time dilation fac¬ 
tor is calculated according to 7 = (1 — In order 

to reach sufficient accuracy, the calibration reported in 


[13] is applied for the ISOLDE high voltage and the post¬ 
acceleration voltage is measured repeatedly with a Julie 
Research 10-kV high-voltage divider providing a relative 
accuracy of 1.2 • lO”"*. 


3 Analysis and theory 

Sufficient production rates of > 10® ions/s to perform laser 
spectroscopy measurements on atomic Cd with continuous 
ion beams were reached for i06-i24,i26(^^^ addition to 
the nuclear ground states the I = 11 / 2 “ isomers were ob¬ 
served for all odd-mass isotopes from ^^^Cd upwards. The 
isotope with the highest natural abundance in the cad¬ 
mium chain, ^^'^Cd, was chosen as the reference for isotope 
shift measurements. The spectra of the other isotopes were 
recorded in alternation with ^^^Cd. A typical spectrum of 
the reference isotope is shown in Fig. 1 as a function of 
the post-acceleration voltage at the charge-exchange cell. 
It exhibits an asymmetric tail at lower voltages, i.e. higher 
beam energies. Thus it is produced by ions in the beam 
that need additional post-acceleration to be in resonance 
because they suffer a loss of energy in atomic collisions be¬ 
fore they interact with the laser beam. The spectrum can 
be well fitted by two Voigt profiles with identical line shape 
parameters. The complete set of reference spectra were 
first analyzed using these Voigt profiles with free distances 
and relative intensities. It was observed that the distance 
between the two peaks was approximately constant with 
an average value of AU = 13.5V but the common width 
and relative intensity changed slightly during the four days 
of beam time. Such a tail is often observed in collinear 
laser spectroscopy and usually ascribed to energy loss in 
the charge-exchange reaction or ion-atom collisions along 
the beam line. However, in this case the uncommonly large 
separation is not explainable by a single process. Instead, 
multiple side peaks shifted in the energy by integer mul¬ 
tiples of the transition energy would give a similarly good 
description of the line shape. Hence, the single side peak 
at this large distance can probably be understood as a 
feature effectively mimicking the resulting tail from mul¬ 
tiple collisional processes. To avoid artificial shifts of the 
line center in the fitting routine, all spectra - including 
those of the reference isotope - were finally analyzed with 
the double Voigt profile of common linewidth and a fixed 
peak distance. The reference isotope was fitted with free 
intensity ratios. The spectra of all other isotopes were an¬ 
alyzed twice: once with linewidth parameters and relative 
intensities as obtained from the preceding reference scan 
and once with parameters as determined in the succeed¬ 
ing reference scan. The resonance positions and hyperfine 
parameters obtained in the two fits were averaged and the 
larger of the two fitting uncertainties was adopted. 

The recorded hyperfine spectra of the odd isotopes 
were fitted using the centre of gravity (eg) and the A 
and B factors of both electronic levels as free parame¬ 
ters. The hfs shift of each sublevel with total angular mo¬ 
mentum F = I + J composed of electronic angular mo¬ 
mentum J and nuclear spin I and corresponding factor 
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Fig. 1. Spectrum of the reference isotope ^^^Cd in the 
5s 5p ®P 2 —>■ 5s 6s ®Si transition. The fitted curve is the sum of 
two Voigt profiles - also shown individually - with a relative 
position of 13.5 V to each other. The lower panel shows the 
residuum between the experimental data and the fit. Typical 
total linewidth (Lorentzian + Gaussian contribution) is about 
30-40 MHz. 


C = F{F -b 1) — /(/ -I- 1) — J{J + 1) was calculated ac¬ 
cording to the first order hyperfine energy 


= A— -b B 


C(C' + l)-/(/ + l)J(J + l) 


2/(2/- 1)J(2J- 1) 


■ ( 8 ) 


In the fit function, these shifts determine the spectral line 
positions relative to the centre of gravity and each hyper¬ 
fine component is again composed of the main peak and 
the tailing peak as discussed above. 

The hyperfine parameters A and B of the upper and the 
lower electronic levels were fitted independently for all iso¬ 
topes except for ^^®Cd for which the spectrum is shown 
in Fig. 2. In this isotope, due to overlapping resonances 
of the isomer and the nuclear ground state, shown and 
enlarged in Fig. 2, the ratio of the two A factors was fixed 
to the error-weighted mean value of the other isotopes as 
discussed in Sect. 4.2. However, in an independent anal¬ 
ysis, the spectrum could also be analyzed without this 
constraint by fitting both nuclear states simultaneously 
and both results were in agreement within uncertainties. 

The uncertainty of the acceleration voltage gives rise to 
systematic uncertainties of the determined isotope shifts. 
Compared to these, the contributions from atomic masses 
and the laser frequency are negligible. 


4 Results and discussion 

The fitting procedure yields the hyperfine parameters A 
and B and the centroid transition frequencies Vcg, relative 
to ^^^Cd. From these parameters, the nuclear observables 
can be extracted if the hyperfine fields are known at the 
site of the nucleus. This section describes how these atomic 
parameters were determined. 


4.1 Nuclear spin 

The nuclear spin of an isotope can be directly deduced 
from the number of peaks if / < J (in the case I = 1/2 also 
for / = J). In the 5s5p^P2 (J = 2) —)• 5s6s ^Si(J' = 1) 
transition, for example, this determination of the nuclear 
spin is possible for / = 1/2 and / = 3 / 2 . In these cases three 
peaks and eight peaks are observed, respectively (the case 
1 = 1 cannot occur for the Cd isotopes). For larger spins 
nine hyperfine components appear. The situation is illus¬ 
trated for the ground state of ^^®Cd (dark blue) and its 
11 / 2 “ isomer (light blue) in Fig. 2. Both these states can 
be clearly distinguished from each other. Nine peaks be¬ 
long to the 11 / 2 “ isomer, exhibiting the same hyperfine 
pattern as the other isomers. Only three peaks remain for 
the ground state of ^^®Cd which therefore must have a 
nuclear spin of / = 1 / 2 . This clarifies a misassignment in 
literature where ^^®Cd was assigned a spin of / = 3/2 based 
on 7 -ray studies [14]. This is consistent with conclusions 
drawn from the measurements on Cd'*" ions published pre¬ 
viously [10]. For all other isotopes, the formerly (partly 
tentatively) reported nuclear spins were confirmed by the 
hfs measurements. The nuclear spins are listed in Table 3. 
Revisions of spin assignment occur at ^^^Cd from I = 5/2 
to / = 1/2 and at ^^^Cd from / = 1/2 to / = 3/2 for the odd 
isotopes. All isomers have the spin I = 11 / 2 . According to 
the spins, the unpaired neutron should occupy the 2d5/2, 
3 si/ 2 , 2d3/2 and I/ 111/2 shell-model orbitals, respectively. 


4.2 Hyperfine structure parameters 

The extracted hfs parameters A and B of the studied 
transition are presented in Table 1. The quoted uncertain¬ 
ties include the statistical and the systematic uncertainties 
added in quadrature. The latter are caused by the uncer¬ 
tainty in the determination of the post-acceleration volt¬ 
age applied to the charge-exchange cell. The quadrupole 
interaction for a pure 5s6s^Si state should be zero. How¬ 
ever, since the total angular momentum is J = 1, con¬ 
figuration mixing with other J = 1 states results in a 
(usually small) finite B value. Therefore, the B param¬ 
eter for the 5s6s^Si electronic level was not fixed to be 
zero in the fitting routine. The resulting B factors of this 
level are nevertheless very small, have uncertainties much 
larger than their magnitudes and are thus not included in 
Table 1. 

In Fig. 3 the ratio of the A factors A{5s5p ^P 2 ) and 
A(5s 6s ^Si) is plotted. The values of the ground and iso¬ 
meric state of ^^®Cd are missing since in these cases the 
ratio was constrained in the fitting as discussed above. The 
ratio is constant within uncertainties, showing no appre¬ 
ciable amount of isotope or spin dependence - a fact that 
supports the choice of an average A-factor ratio for fitting 
the spectrum of ii9Ti9™Cd. The error-weighted mean A 
factor ratio of all isotopes is 0.42479(3), indicated by a 
dashed line in Fig. 3. 

Other results available from literature are included in 
Table 1. They are all in agreement with the values ob¬ 
tained in this work. While A factors for the 5s 5p ^P 2 
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Fig. 2. (color online) Hyperfine structure of the nuclear ground state (dark blue) and the isomeric ii/ 2 “state (light blue) of 
^^®Cd in the 5s 5p ®P 2 —>■ 5s 6s ®Si transition. Both states can be clearly distinguished from each other. The number of hyperfine 
lines arising from the nuclear ground state unambiguously yields a nuclear spin of / = 

Table 1. A and B factors of the hfs splitting of the 5s 5p ®P 2 and the 5s 6s ®Si electronic levels in the cadmium atom for 
isotopes with mass number A = Z -\- N, half-life T 1/2 and spin and parity 7'^. The B factor of the upper electronic level is within 
its uncertainty compatible with zero and is not included in the table. The listed total uncertainties are the quadratic sum of the 
statistical uncertainty and the systematic uncertainty as discussed in the text. For ^^®Cd and ^^®™Cd the ratio of the A factors 
was constrained to the average of all other isotopes. Literature values are listed where available with their references (square 
brackets) for comparison. 





5s 5p ®P 2 

5s6s ®Si 

5s 5p ®P 2 

AlNMR [18] 

Z + N 

7i/2 

r 

A (MHz) 

A (MHz) 

B (MHz) 

(m) 

107 

6.50 h 

5 / 2 + 

-682.0 (4) 

-1605.6 (7) 

-7290 (3) 

-0.6150554(11) 

109 

461.4 d 

^l2+ 

-917.5 (1) 

-2159.9 (4) 

-7296 (2) 

-0.8278461(15) 

111 

stable 

V 2 + 

-3292.7 (5) 

-7750.9(10) 


-0.5948861 (8) 




-3292.9364 (2)® 

-7750.5(18^ 



113 

stable 

Y 2 + 

-3444.5 (6) 
-3444.6344(16^ 

-8108.8(11) 


-0.6223009 (9) 

115 

53.46 h 

Y 2 + 

-3588.6 (8) 

-8444.3(28) 


-0.6484259(12) 

117 

2.49 h 

Y 2 + 

-4116.7 (8) 

-9687.9(17) 






-4116.1(63Y 

-9698(10Y 



119 

2.69 min 

V 2 + 

-5091.8(12) 




121 

13.5s 

3 / 2 + 

-71156.4 (4) 

-72722.0 (6) 

-141 (7) 


123 

2.1s 

3 / 2 + 

-71457.2 (8) 

-73431.9(13) 

+31 (4) 


111m 

48.50 min 

IV 2 " 

-556.5 (1) 

-1310.3 (2) 

-363 (2) 

-1.1051 (4) 

113m 

14.1a 

IV 2 " 

-547.8 (1) 

-1289.5 (2) 

-299 (2) 

-1.087784 (2) 

115m 

44.56 d 

IV 2 " 

-524.3 (2) 

-1234.2 (2) 

-234 (2) 

-1.0410343(15) 

117m 

3.36 h 

IV 2 " 

-502.0 (1) 

-1181.8 (3) 

-156 (1) 





-501.4 (6Y 

-1181.5 (8Y 

-144(15Y 


119m 

2.20 min 

172 " 

-485.3 (2) 


-67 (1) 


121m 

8.3 s 

IV 2 " 

-508.5 (2) 

-1197.2 (3) 

+3 (3) 


123m 

1.82s 

IV 2 " 

-504.0 (3) 

-1186.5 (5) 

+63 (6) 



^ Faust et al. [15] 

^ Brimicombe et al. [16] 
" Boos [17] 
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Fig. 3. Ratio of the A factors of the 5s 5p ^P 2 and the 
5s 6s ^Si electronic levels obtained in the nuclear ground state 
(•) and the isomeric (a) state. The dashed line indicates the 
error weighted mean value and the dotted lines represent the 
la uncertainty. 

level in and ^^^Cd (both I = 1 / 2 ) are available in 

the literature with orders of magnitude higher precision 
[15], obtained with the atomic beam magnetic resonance 
technique [15], corresponding values for the other isotopes 
are from pressure-scanned Fabry-Perot spectroscopy [16] 
and from collinear laser spectroscopy performed at COL- 
LAPS in 1989 [17]. Both have larger or comparable un¬ 
certainties than the values reported here. Accurate NMR 
values for the nuclear magnetic moment are available for 
several isotopes and are included in Table 1 since they are 
used in the following to analyse hfs anomalies and serve 
as a reference for the extraction of nuclear moments. 


4.3 Hyperfine structure anomaly and the 
Moskowitz-Lombardi rule 

A change of the magnetization distribution from one iso¬ 
tope to the other gives rise to an (isotopic) hfs anomaly 

1A2 

as defined in Eq. (5). For the isotopes with precisely 
known nuclear magnetic moments obtained by NMR, the 
(isotopic) hfs anomalies for the isotopes with mass 

number^ M with respect to ^^^Cd are calculated from the 
fitted hfs parameters of the initial and the final electronic 
level A(^P 2 ) and A(^Si), respectively. These anomalies 
are listed in Table 2 and are compared to values from the 
literature for the 5s 5p ^Pi level. 

Even though the ratio of the A factors in the upper and 
the lower level of the transtion as plotted in Eig. 3 is con¬ 
stant within the measurement uncertainty, the individual 
levels exhibit a substantial amount of hfs anomaly. This 
is however masked by the fact that the so-called differen¬ 
tial hfs anomaly, being the difference in the hfs anomaly 

^ Here we use M instead of the usually used letter A for 
the mass number of an isotope to avoid confusion with the 
hyperfine A factor 


between the initial and the final level of the atomic tran¬ 
sition, is very small. This leads to the conclusion that the 
hfs anomaly arises dominantly from the 5s electron which 
is common to both atomic levels. The 5p electron in the 
^P 2 level must have largely 5 p 3/2 character to build up 
the J — 2 level and thus has a probability density van¬ 
ishing in the nuclear interior. Consequently, also the 6 s 
electron cannot have a considerable contribution, other¬ 
wise the differential hfs anomaly between the two elec¬ 
tronic levels would be larger. A comparison between the 
observed isotopic hfs anomaly in the 5s 5p ^P 2 level and 
the one in the 5s 5p ^Pi fine structure level (lower rows of 
Table 2) indicates that the hfs anomaly in the ^Pi level is 
even slightly smaller. Since this fine structure level must 
have some p 1/2 character, the p 1/2 contribution must be 
responsible for the difference and seems to compensate 
partially the hfs anomaly caused by the 5s 1/2 electron. 

This line of reasoning can be supported by relating 
the A factors in the two-electron systems of the cadmium 
atom to the single-electron a-factor in the ion. This can be 
performed assuming Russel-Saunders (LS) coupling and 
neglecting contributions from the 5p electron. This leads 
to [19] 


A (5s 6s ^Si) Rs i (oss -1- aes) 

(9) 

r 3r) \ ^ 

5P P2)“»=-2(f+l)= 4 ^ 

( 10 ) 


Using the A and g factors of the ion [10] and the atomic 
levels observed here, we find 

A/g (5s ^ 81 / 2 ) ~ 12, 200 MHz, (11) 

A/g (5s 6 s ^Si) Ri 6 , 500 MHz, and (12) 

A/g {5s 5p ^¥ 2 ) ~ 2,800 MHz. (13) 

Based on the first value and the above quoted relation 
(10) one would expect an A/g-value for the 5s 5p ^P 2 
level of « 3,050 MHz which is even 10% larger than the 
experimental value, indicating that indeed the hyperfine 
interaction arises almost exclusively from the 5s electron, 
reduced by a weak shielding from the 5p electron. Using 
Eq. (9), it follows that the 6 s electron has a contribution 
that is of the order of 10% of the 5s electron. 

The isotopic hfs anomaly for the isotope pairs 
and is close to zero in all electronic levels. This 

is expected since these isotopes have identical nuclear spin 
1 = 1/2 and should therefore have a similar magnetization 
distribution. Eor the isotopes with a larger nuclear spin, 
larger values are observed but again isotopes having iden¬ 
tical nuclear spin and parity, e.g. io'i',i09(^j (V^"*") 01 the 
11 / 2 “ isomers, exhibit niA^ values that are similar in 
size. 

Since we have values for the hfs anomaly of several 
different nuclear states, it is worthwhile to check whether 
the Moskowitz-Lombardi rule that has been established in 
the lead region [2, 3], works also well in the tin region. Eor 
Cd, all isotopes for which the hyperfine anomaly can be 
determined are associated with si/ 2 , £^ 5/2 and hii /2 shell- 
model states and have nuclear spin / = £ -|- 1 / 2 . To ex¬ 
amine whether a similar linear behavior holds for the Cd 
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Table 2. Isotopic hyperfine anomalies m of isotopes with 
mass number M and spin I relative to ^^^Cd in the electronic 
levels 5s 5p ®P 2 and 5s 6s ®Si. Listed are the isotopes for which 
high-precision values of the nuclear magnetic moment are avail¬ 
able from literature. The hyperfine parameters A obtained in 
this work, which were used for the calculations according to 
Eq. (5), are listed in Table 1. Literature values for the 5s 5p ®Pi 
level are included for discussion. 


Isotope 


5s 5p 3p2 

5s6s 


M 

r 

(%) 

“iA*f [%] 

Ref. 

107 

5/2+ 

-0.17(5) 

-0.18(4) 

this work 

109 


-0.12(1) 

-0.12(1) 

this work 

113 

1/2+ 

0 .00(1) 

-0.01(1) 

this work 

113 

1/2+ 

-0.00143(6) 

-0.01(4) 

[15, 20] 

115 

1/2+ 

0 .01(2) 

0.05(3) 

this work 

111 m 

11 / 2 " 

-0.08(4) 

-0.10(4) 

this work 

113m 

11 / 2 " 

-0.08(2) 

-0.08(1) 

this work 

115m 

11/2" 

-0.09(4) 

-0.09(2) 

this work 

Isotope 


5s 5p ®Pi 



M 

r 

(%) 


Ref. 

107 

5/2+ 

-0.0958(8) 


[20, 21] 

109 

5/2+ 

-0.0912(7) 


[20, 21] 

113 

1/2+ 

-0.00023(40) 


[7, 20] 

113m 

11 / 2 " 

-0.0773(5) 


[7, 20] 



Fig. 4. Moskowitz-Lombardi rule for Cd isotopes: the hfs 
anomaly for the 5s 6s ®Si level is plotted as a function of the 
difference of the inverted nuclear moments for nuclear states 
with I = £-1- 1 / 2 . A common linear relation is only observed for 
states with 1^0. 


isotopes, the isotopic hfs anomaly is plotted in Fig. 4 as 
a function of the difference of the inverted nuclear mo¬ 
ments for the case of the 5s 6s ^Si level. The distribution 
looks very similar for 5s 5p ^P 2 and it is obvious that 
those isotopes that have a ds/a or hiy^ character do in¬ 
deed exhibit a common linear relation, whereas the iso¬ 
topes with an si /2 configuration show a distinctly differ¬ 


ent behavior. The s states belong neither to the i + 1/2 
nor to the £ — 1/2 branch, because there is no angular 
momentum £ to which the spin could couple. Whether 
the different isotopes with this ground state configuration 
follow a similar rule with a different slope cannot be de¬ 
cided with the available data. In the case of Hg, there 
is no similar case since no Si /2 isotopes exist in the shell 
with negative parity filled between N = 82 and N = 126. 
The linear fit for the (£ 5/2 and hii /2 states of Cd produces 
a downsloping line with a = —0.135(16) • pN- This 
slope is almost an order of magnitude smaller and of op¬ 
posite sign than in the case of Hg. The corresponding fit 
for the anomalies extracted from the 5s 5p ^P 2 level yields 
a slope of —0.124(33) • pn that agrees well with the 
one of the 5s 6s ^Si level - as expected since the differen¬ 
tial hfs anomaly almost vanishes. Analyzing the available 
data from literature in the 5s 5p ^Pi level, a much smaller 
slope of only a = —0.027(8) • 10~‘^ pN is obtained. Finally, 
our ionic data presented in [10] for the 5s ^Si /2 level result 
in a slope of a = —0.065(38) • 10“^ pN, being roughly 50% 
of the slope for the atomic levels investigated in this work. 
In all cases the isotopes having nuclear s states must be 
excluded to observe a linear behavior. It should be noted 
that the y-axis intercept is in all cases compatible with 
zero within the corresponding uncertainty. 

4.4 Nuclear magnetic dipole moments 

To determine the nuclear magnetic moments we use Eq. (1). 
The hfs anomaly is negligible if the reference dipole mo¬ 
ment is taken from an isotope with the same nuclear spin / 
parity state, as it has been discussed in the previous sec¬ 
tion. This is possible for 5 / 2 +, 1 / 2 + and 11 / 2 “ isotopes, 
where we use the reference isotopes ^°®Cd, ^^^Cd and 
respectively. The corresponding magnetic moments 
from literature that were used are included in Table 3 
and marked with For the two isotopes with nuclear 
spin / = 3 / 2 +, ^^^Cd and ^^^Cd, the magnetic moment 
of ^^^Cd was used as a reference. To estimate the effect 
of a potential hfs anomaly, the largest value of the de¬ 
termined anomalies (0.17%) in Table 2 was assumed and 
the corresponding effect was added quadratically to the 
statistical uncertainty of the respective isotope. The mag¬ 
netic moments listed for Cd I in Table 3 are the weighted 
mean values obtained from the A factors of the ^P 2 and 
the ^Si levels, being consistent within uncertainty. They 
agree with the moments measured by NMR and confirm 
the measurements on ionic cadmium in [10]. 

4.5 Electric field gradient of the 5s5p ^P 2 level in Cd I 

The quadrupole term of the hfs formula (8) contains the 
parameter B which describes the interaction of the nuclear 
electric quadrupole moment Qgwith the electric field gra¬ 
dient 142 produced by the electron cloud at the site of the 
nucleus, according to 


B = eQs Fzz 


(14) 
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Table 3. Nuclear magnetic dipole moments and electric qnadrupole moments of the cadmium isotopes. For the calculation of the 
magnetic moments different reference values (marked with are used depending on the nuclear spin. The magnetic moments 
listed for Cd I are the weighted mean values obtained from the A factors of the ®P 2 and the ®Si levels. The columns labeled 
Cd'^" show the magnetic moments and the quadrupole moments obtained from collinear laser spectroscopy of ionic cadmium in 
the 5s ^Si /2 —>■ 5p ^P 3/2 transition reported in [10]. They are in reasonable agreement with the atomic data. The reference value 
for the electric quadrupole moments is the qnadrupole moment of ^°®Cd from [10] (reasons for that are explained in Sect. 4.6). 
The second parentheses for the Q-values indicate the systematic uncertainty from the uncertainty of the electric field gradient. 





Cd I 

Cd+ [10] 

5s5p ^P2 

Cd+ 5p "P 3/2 [10] 

A 

r 

AtNMR (Pn) 

p (Pn) 

P (Pn) 

Qs (b) 

Qs (b) 

107 

^12+ 

-0.6150554(11) 

-0.6154 (3) 

-0.6151(2) 

+0.593 (7)(25) 

+0.601 (3) (24) 

109 

5/2+ 

-0.8278461(15)# 




+0.604(1)(25) # 

111 

1 / 2 + 

-0.5948861 (8)# 





113 

1 / 2 + 

-0.6223009 (9) 

-0.6224 (1) 

-0.6224(2) 



115 

1 / 2 + 

-0.6484259(12) 

-0.6483 (2) 

-0.6483(2) 



117 

1 / 2 + 


-0.7437 (1) 

-0.7436(2) 



119 

1 / 2 + 


-0.9199 (2) 

-0.9201(2) 



121 

3 / 2 + 


-t0.6268(ll) 

-t0.6269(7) 

-0.288(15)(12) 

-0.274(7)(11) 

123 

3 / 2 + 


-t0.7900(14) 

+0.7896(6) 

+0.063 (7) (3) 

+0.042(5)( 2) 

111m 

11/2" 

-1.1052 (2) 

-1.1052 (2) 

-1.1052(3) 

-0.742 (5)(31) 

-0.747(4) (30) 

113m 

11/2" 

-1.0877842(17) 

-1.0877 (2) 

-1.0883(3) 

-0.612 (4)(25) 

-0.612(3)(25) 

115m 

11 / 2 " 

-1.0410343(15)# 



-0.477 (4)(20) 

-0.476(5)(19) 

117m 

11/2" 


-0.9969 (2) 

-0.9975(4) 

-0.319 (2)(13) 

-0.320(6)(13) 

119m 

11 / 2 " 


-0.9636 (5) 

-0.9642(3) 

-0.136 (3) (6) 

-0.135(3)( 5) 

121m 

11/2" 


-1.0098 (3) 

-1.0100(4) 

+0.007 (6) (1) 

+0.009(6) 

123m 

11/2" 


-1.0008 (5) 

-1.0015(3) 

+0.128(11) (5) 

+0.135(4)( 6) 


This electric field gradient (EFG) has been evaluated 
for the 5s5p ^P 2 excited level of neutral cadmium within 
the multi-configuration Dirac-Hartree-Fock (MCDHF) the¬ 
ory [22]. Numerical-grid wave functions have been gener¬ 
ated by means of the atomic structure code GRASP [23, 
24] as self-consistent solutions of the Dirac-Hartree-Fock 
equations [25]. 


4.5.1 Multi-configuration Dirac-Hartree-Fock theory 

In the MCDHF theory, an atomic level is approximated 
by a linear combination of configuration state functions 
(CSF) <?( 7 fc J) of the same symmetry 

NCF 

'f'a(J) = '^Ck{a)<P{jkJ), (15) 

k 

where NCF refers to the number of CSFs and {ck{a)} to 
the representation of the atomic level in the given basis. 
Moreover, the sets jk describe all quantum numbers that 
are required to distinguish the basis states uniquely. In the 
GRASP code, the CSFs are constructed as anti¬ 

symmetrized products of a common set of orthonormal 
orbitals, and are optimized together on the basis of the 
Dirac-Coulomb Hamiltonian 

Hoc = X! ■ Pi + ^(^i)] + X! • 

i i>j 

(16) 


Further relativistic effects due to the Breit interaction 
could, in principle, be added to the representation {cfc(a)} 
by diagonalizing the Dirac-Coulomb-Breit Hamiltonian 
matrix but they are known to have little effect upon the 
EFG and hyperfine parameters of light- and medium-Z ele¬ 
ments [26-28]. For the 5s5p ^P 2 level of neutral cadmium 
only CSFs of even parity and total angular momentum 
J = 2 need to be taken into account in the wave function 
expansion. 


4.5.2 Wave function generation 

The computational methodology of generating wave func¬ 
tions for hyperfine calculations have been explained else¬ 
where [29,30]. We shall describe here only the basics of 
the models which were applied in the three series of com¬ 
putations. All three computational models share a com¬ 
mon set of spectroscopic (occupied) orbitals. All spectro¬ 
scopic orbitals Is, 2s, ..., Ad, 5s, 5p of the 5s5p ^P 2 level 
of neutral cadmium were obtained in the Dirac-Hartree- 
Fock approximation and were frozen in all further steps of 
the computations. Properly antisymmetrized linear com¬ 
bination of products of the set of spectroscopic orbitals 
formed the zero-order representation (often referred to as 
a reference configuration state function) of the 5s5p ^P 2 
atomic level. 

Then two sets of virtual orbitals were generated in a 
series of steps in order to incorporate electron-electron cor¬ 
relations at increasing level of (computational) complexity 
and to monitor their effect upon the level energy and the 
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EFG of the 5s5p ^P 2 level. Our computational models dif¬ 
fer in the way the virtual orbitals are generated, as well as 
in the sets of CSFs in the wave function expansion (15); 
see Refs. [30,31] for further information. Below, we briefly 
refer to these models as single substitutions (S) and sin¬ 
gle with restricted double substitutions {SrrD and SrD); in 
some more detail, these models include: 

S Single substitutions from the valence and core or¬ 
bitals to eight layers of virtual orbitals. Within 
hrst-order perturbation theory, only single substi¬ 
tutions contribute to the hyperfine energy [32,33]. 
In the MCDF method, in contrast, double substi¬ 
tutions from the (sub-) valence shells often domi¬ 
nate energetically over correlation corrections and, 
hence, the virtual orbitals obtained in correlated 
variational calculations are determined predomi¬ 
nantly through the effects of double substitutions 
upon the total level energy. For this reason, how¬ 
ever, it is important that the virtual space is con¬ 
structed from orbitals optimized for single substi¬ 
tutions [34]. 

SrrD Single and restricted double substitutions from the 
valence and core orbitals to seven layers of virtual 
orbitals. Double substitutions were restricted in 
the sense that at most one electron is promoted 
from core orbitals, while the other (or both) is 
(are) promoted from the valence shells 5s5p. The 
size of the multi-configuration expansions is fur¬ 
ther restricted in this model by eliminating all 
CSF which do not interact with the reference CSF. 
Calculations based on an SrrD computational mo¬ 
del have been found optimal for ab initio calcula¬ 
tions of hyperfine structures in heavy atoms [30]. 

SrD The orbital sets as generated in the step SrrD 
above were utilized also for performing configura¬ 
tion-interaction calculations (i.e. with all orbitals 
frozen), where the CSF reductions described in the 
step SrrD above were lifted. For such expansions, 
the SrD approximation can be considered as a con¬ 
sistency check on the wave functions obtained in 
the SrrD model. 

For the two models S and SrrD, virtual orbitals were gen¬ 
erated stepwise in terms of layers [30,35], until satura¬ 
tion of the calculated electric field gradient (FFC) was 
observed. Fight layers of virtual orbitals were generated 
in the approximation S by opening eventually all core or¬ 
bitals, down to Is, for substitutions. For the SrrD and SrD 
models, a total of seven layers of correlation orbitals were 
generated. 

In ah initio calculations of hfs, double and triple sub¬ 
stitutions are known to cancel each other to a large de¬ 
gree [34], and their combined effect upon the calculated 
EFG values is typically a few percent or even less. This 
was confirmed by a number of earlier calculations on other 
medium and heavy elements [30,31,36]. Such unrestricted 
double and triple substitutions were included also in our 
recent calculations for the excited 5p level of ^°®Cd''" 
ion [10], in which these double and triple substitutions 
contributed individually less than 3 percent each, and be¬ 


low 1 percent when considered together (i.e. well below 
the estimated 4 percent accuracy of the Hnal EFG value 
in Ref. [10]). In the present work, we therefore performed 
only three large-scale configuration interaction calcula¬ 
tions, in which we compared the effects of unrestricted 
double against unrestricted triple substitutions. These com¬ 
putations were carried out for (i) substitutions from 4s4p- 
4d5s5p occupied orbitals to one layer of virtual orbitals; 
and (ii) for substitutions from the 3s3p3d4s4p4d5s5p oc¬ 
cupied orbitals to one layer. The results of these three 
conhguration-interaction calculations are represented by 
the three unconnected (magenta) circles on the right side 
of figure 5. They clearly indicate that the corrections aris¬ 
ing from the unrestricted double and triple substitutions 
are smaller than the estimated 4 percent accuracy of the 
final EFG value. 

The wave functions have been generated with the nucleus 
modelled as a sphere, and a two-parameter Fermi distri¬ 
bution [37] was employed to approximate the radial de¬ 
pendence of the nuclear charge density. All other nuclear 
electromagnetic moments were assumed to be point-like, 
i.e. the magnetization distribution inside the nucleus (the 
Bohr-Weisskopf effect [1]), and still higher nuclear elec¬ 
tromagnetic moments, as well as Breit and QED correc¬ 
tions [38], were all neglected, since they are expected to 
be negligible at the current level of accuracy. These effects 
were found to be small even for an element as heavy as 
radium [39,40]. All further details of the computations, in¬ 
cluding a method to take advantage of the computational 
balance between the hyperfine A and B coefficients in the 
evaluation of EFG, will be described elsewhere. 

4.5.3 EFG theory 

The quadrupole moment Q of an isotope is related to the 
(electric-quadrupole) hyperfine constant B through the 
equation 

B(J) =2eQ , (17) 

with e being the charge of the electron, and where the 
(second-rank) operator acts upon the electronic co¬ 
ordinates of the wave functions. It is this operator in the 
theory of hfs [32,33] that describes the electric field gra¬ 
dient, EFG, at the site of the nucleus 

(jj|r<qjj) = i (0) , (18) 

4.5.4 EFG results 

Figure 5 displays the calculated values of the EFG for 
the 5s5p ^P 2 excited level of neutral cadmium for dif¬ 
ferent computational models at various degrees of com¬ 
plexity. Results are shown for the three (computational) 
models as explained above, i.e. by including: single substi¬ 
tutions S (blue diamonds), single with restricted double 
substitutions SrrD (olive squares), as well as for the SrD 
conhguration-interaction calculations (magenta circles). 
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number of virtual orbital layers in MCDHF expansion 



number of virtual orbital layers in MCDHF expansion 

Fig. 5. (Color online) Electric field gradient eVzz (top) and 
magnetic field {A/g, bottom) for the low-lying excited 5s5p ®P 2 
level of neutral cadmium and for the three computational mod¬ 
els as described in the text, i.e. by including single substitutions 
S (blue diamonds); single with restricted double substitutions 
SttD (olive squares); and StD expansions (magenta circles). 
The integers on the abscissa axis represent the number of (cor¬ 
relation) layers of virtual orbitals that are taken into account 
in the expansion (15). The three (magenta) circles to the right 
side of the figure display configuration-interaction calculations 
which include full double and triple substitutions. All lines are 
drawn only for the guidance of the eyes. The two (red) horizon¬ 
tal lines without symbols in the top graph denote the estimated 
uncertainty of the final computational result of el4z, which is 
indicated by the dashed line. In the bottom graph, the (red) 
horizontal line represents the experimental value of A/g. See 
text for further details. 


All lines in this graph are drawn only for the guidance 
of the eyes. The computational complexity is represented 
by the number of (correlation) layers on the abscissa axis 
that are taken into account in the expansion (15) of the 
corresponding wave functions. Up to 3868 CSFs were in¬ 
cluded in the expansions for single substitutions into eight 
correlations layers, and up to 216643 CSFs in the StD 
model in the final stages of the computations. The uncor¬ 
related Dirac-Hartree-Fock value (n = 0) is also shown 


but is not considered for estimating the uncertainty of 
the computations. Our largest expansions yield the value 
eVzz = 500 ± 30 MHz/b, where the uncertainty was as¬ 
sumed to cover all values for the layers of virtual orbitals 
with n = 2,...,8. This conservative estimate of the ‘un¬ 
certainty’ due to missing correlations is shown by the two 
horizontal lines (without symbols) in Fig. 5. 

While the final value el4z = 500.108 MHz/b from the 
largest computation almost coincides (perhaps acciden¬ 
tally) with the n = 7 points of the SrrD and SrD ap¬ 
proaches in Fig. 5, single substitution model S yields the 
final value close to the upper straight (uncertainty) line. 
As mentioned above, the two models SrrD and SrD are 
considered to be suitable especially for ab initio calcula¬ 
tions of hyperfine structures. Therefore, it appears well 
justified to assume the last points of the SrrD and SrD 
curves as final value, and the deviations of all three curves 
together as a measure for the uncertainty of the compu¬ 
tations. 

As a test of their reliability, the same wave functions 
were used to calculate the magnetic field produced by the 
electrons at the site of the nucleus, and subsequently to 
compare with the measured A-factors and magnetic mo¬ 
ments. The bottom graph in Fig. 5 displays the calculated 
values of the magnetic field A/g for the 5s5p ^P 2 excited 
level of neutral cadmium, obtained from the same compu¬ 
tational models as described in Sec. 4.5.2 above and in the 
caption of Fig. 5. The horizontal (red online) straight line 
represents the experimental value of A/g. These graphs 
support the conclusions drawn in the preceeding para¬ 
graphs, concerning the uncertainty of the computations. 
The graphs clearly indicate that the calculations of the 
magnetic field have not produced a converged value of 
A/g, but A/g is not the objective of this work. 

From the above analyses, we conclude a (theoretical) 
value eVzz = 500 ± 30 MHz/b, which can be utilized in 
deriving the quadrupole moments from the measurements 
of the hyperfine B parameters. 

4.6 Electric quadrupole moments 

The nuclear quadrupole moments can now be directly ex¬ 
tracted from the determined spectroscopic B factors using 
the calculated EFG value of the 5s5p ^P 2 level. For ^°®Cd, 
we obtain Qs(^°®Cd) = 592(l)(36)b. The first paranthe- 
ses represent the statistical measuring uncertainty and the 
second ones the EFG-related systematic uncertainty. The 
latter is about 50% larger than in the value (5s(^°®Cd) = 
0.604(1) (25) b, reported in [10] from the 5p^P3/2 level of 
the Gd’*' ion and based on eVzz = 666(27) MHz/b. This 
EFG for the Gd'*' level was calculated in the same way 
as described in the previous section. The quadrupole mo¬ 
ment extracted directly from the atomic level is only 2% 
smaller than the result obtained for the ion, a difference 
that is insignificant at the level of accuracy of the EFG 
values. Hence, the results obtained in two different ways 
clearly demonstrate that the EFG calculations for the 
atomic and the ionic systems are compatible and pro¬ 
vide consistent data. In order to provide a coherent set 
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of quadrupole moments from both measurements with the 
smallest possible uncertainty, we prefer to base all spectro¬ 
scopic quadrupole moments on the ionic result by using 
Qs(^°®Cd) = 0.604(1)(25)b as reference value and em¬ 
ploying the relation 

for the other isotopes. The results are listed in Table 3. 
The most intriguing observation is the linear increase of 
the electric quadrupole moments of the 11 / 2 “-isomers that 
has been extensively discussed in [10]. The values reported 
here are in excellent agreement with the ionic results. It 
should be noted that using the directly calculated EFG 
value for the atomic level would simply result in a reduc¬ 
tion of all quadrupole moments by 2%. 

The literature values [18] of the quadrupole moments 
of lO^Cd, l°^Cd, l°®Cd, l“™Cd, 

more than 10 % larger than the present results. They are 
all based on the hfs of the 5s 5p ^Pi level of Cd I and 
a semi-empirical evaluation of the EFG using information 
from the magnetic hyperfine interaction and the fine struc¬ 
ture splitting [9] and neglecting Sternheimer-type shield¬ 
ing correction. A systematic uncertainty from this proce¬ 
dure was estimated to be 10 %, which is in accordance 
with the deviation from our results. 

Finally, the isotope shifts with respect to ^^'^Gd and 
the isomer shifts between the isomeric and the respective 
ground state are extracted from the centers of gravity ob¬ 
tained from the hfs fits. These will be discussed in a sep¬ 
arate publication combined with similar information from 
the ionic 5s^Si/2 —>■ 5p^P3/2 transition. 

5 Summary 

Nuclear spins and hfs parameters of and the 

isomers were determined by collinear laser 

spectroscopy. The nuclear spin of ^^®Gd was unambigiously 
determined as / = 1 / 2 . Nuclear magnetic dipole moments 
and electric quadrupole moments were extracted using ref¬ 
erence values and results of EFG calculations, respectively. 
Hyperhne structure anomalies and differential hfs anoma¬ 
lies were investigated and the validity of the Moskowitz- 
Lombardi rule was discussed. While the anomalies of nu¬ 
clear d5/2 and hn/^ states show a linear behavior compara¬ 
ble to that observed in the Hg isotopes, the isotopes with 
an S 1/2 ground state having £ = 0 clearly deviate from the 
rule. The magnetic moments and quadrupole moments are 
in excellent agreement with values measured in the ionic 
system [10]. 


6 Acknowledgement 


This work has been supported by the Max-Planck Society, the 
German Federal Ministry for Education and Research under 


Contract No. 05P12RDCIC, the Helmholtz International Cen¬ 
ter for FAIR (HIC for FAIR) within the LOEWE program 
by the State of Hesse, the Belgian lAP Project No. P7/12, the 
FWO-Vlaanderen, and the European Union seventh framework 
through ENSAR under Contract No. 262010. N. F. received 
support through GRK Symmetry Breaking (DFG/GRK 1581). 
We thank the ISOLDE technical group for their professional 
assistance. The theory part of this work was supported by the 
Polish Ministry of Science and Higher Education (MNiSW) in 
the framework of the project No. N N202 014140 awarded for 
the years 2011-2014. The large-scale calculations were carried 
out with the supercomputer Deszno purchased thanks to the fi¬ 
nancial support of the European Regional Development Fund 
in the framework of the Polish Innovation Economy Opera¬ 
tional Program (contract no. POIG.02.01.00-12-023/08). 


References 

1 . A. Bohr and V. F. Weisskopf, Phys. Rev. 77 , 94 (1950). 

2. P. A. Moskowitz and M. Lombardi, Phys. Lett. B 46, 334 
(1973). 

3. P. A. Moskowitz, Phys. Lett. B 118 , 29 (1982). 

4. J. R. Persson, Hyperf. Int. 162 , 139 (2005). 

5. C. Ekstrom et al., Nucl. Phys. A 348 , 25 (1980). 

6 . A. C. Mueller et al., Nucl. Phys. A 403, 234 (1983). 

7. R. L. Chaney and M. N. McDermott, Phys. Lett. A 29 , 
103 (1969). 

8 . P. W. Spence and M. N. McDermott, Phys. Lett. A 42 , 
273 (1972). 

9. N. S. Laulainen and M. N. McDermott, Phys. Rev. 177 , 
1615 (1969). 

10. D. T. Yordanov et al., Phys. Rev. Lett. 110 , 192501 (2013). 

11. K. Kreim et al., Phys. Lett. B 371 , 97 (2014). 

12. K. Burns and K. B. Adams, Energy Levels and Wave¬ 
lengths of Natural Cadmium and of Cadmium-114, J. Opt. 
Soc. Am. 46 , 94 (1956). 

13. A. Krieger et al., Nucl. Instrum, and Meth. A 632 , 23-31 

( 2011 ). 

14. D. M. Symochko et al., Nucl. Data Sheets 110 , 2945-3105 
(2009). 

15. W. Faust, M. N. McDermott and W. Lichten, Phys. Rev. 
120 , 469 (1960). 

16. M. S. W. M. Brimicombe et al., Proc. R. Soc. Lond. A. 
352 , 141-152 (1976). 

17. N. Boos, Kollineare Laserspektroskopie an neutronen- 
reichen Cadmiumisotopen, Diploma Thesis, Institut fiir 
Kernphysik, Universitat Mainz, 1989. 

18. N. J. Stone, Table of nuclear magnetic dipole and electric 
quadrupole moments. Atomic Data Nucl. Data Tab. 90 , 
75 (2005). 

19. H. Kopfermann, Nuelear Moments, Academic Press, New 
York, 1958. 

20. J. R. Persson, arXiv:1110.5991 (2011). 

21. P. Thaddeus and M. N. McDermott, Phys. Rev. 132, 1186 
(1963). 

22. 1. P. Grant, Relativistic Quantum Theory of Atoms and 
Molecules: Theory and Computation, Springer, New York, 
2007. 

23. P. Jonsson, X. He, and C. Froese Fischer, Comput. Phys. 
Commun. 176 , 597 (2007). 



12 


Nadja Frommgen et al.: Collinear laser spectroscopy of atomic cadmium 


24. P. Jonsson, G. Gaigalas, J. Bieron, G. Froese Fischer, and 
I. P. Grant, Comput. Phys. Commun. 184, 2197-2203 
(2013). 

25. I. P. Grant, Comput. Phys. Commun. 84, 59 (1994). 

26. K. Paduch and J. Bieron, J. Phys. B: At. Mol. Opt. Phys. 
33, 303-311 (2000). 

27. A. Derevianko, Phys. Rev. A 65, 012106 (2001). 

28. J. Bierori, P. Pyykko, D. Sundholm, V. Kello, and A. J. 
Sadlej, Phys. Rev. A 64, 052507 (2001). 

29. J. Bieron, P. Jonsson, and C. Froese Fischer, Phys. Rev. 
A 53, 2181 (1996). 

30. J. Bieron, C. Froese Fischer, P. Indelicate, P. Jonsson, and 
P. Pyykko, Phys. Rev. A 79, 052502 (2009). 

31. J. Bieroii, I. P. Grant, and C. Froese Fischer, Phys. Rev. 
A 56, 316 (1997). 

32. I. Lindgren, Rep. Prog. Phys. 47, 345 (1984). 

33. B. G. Wybourne, Spectroscopic Properties of Rare Earths, 
Wiley, New York, 1965. 

34. B. Engels, Theor. Chim. Acta 86, 429 (1993). 

35. S. Fritzsche, Phys. Scr. TlOO, 37 (2002). 

36. J. Bierori, P. Pyykko, and P. Jonsson, Phys. Rev. A 71, 
012502 (2005). 

37. K. G. Dyall, I. P. Grant, C. T. Johnson, F. A. Parpia, and 
E. P. Plummer, Comput. Phys. Commun. 55, 425 (1989). 

38. B. J. McKenzie, I. P. Grant, and P. H. Norrington, Com¬ 
put. Phys. Commun. 21, 233 (1980). 

39. J. Bierori and P. Pyykko, Phys. Rev. A 71, 032502 (2005). 

40. J. Bierori, J. Phys. B: At. Mol. Opt. Phys. 38, 2221-2228 
(2005). 



